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Electric dipole moment of due to P- and T-violating weak interactions 
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We report rigorous calculations of electric dipole moment (EDM) in ^^®Ra due to parity and 
time-reversal violating tensor-pseudotensor (T-PT) and nuclear Schiff moment (NSM) interactions 
between the electrons and nucleus by employing the relativistic all order coupled-cluster (RCC) 
methods at various levels of approximation. The most accurate EDM (cZa) results are obtained 
as dA = —10.04 X 10“^°Ct < Un > |e|cm and d^^’^ = —6.79 x 10“^^S'(|e|/m®)“'^|e|cm with Ct 
and S are the T-PT coupling constant and NSM respectively. Due to exhaustive treatment of the 
electron correlation effects in these calculations, the EDM results for the corresponding T-PT and 
NSM interactions reduce by about 45% and 23%, respectively, from the previously known values. 
Nonetheless they are still found to be 2-3 times larger than the ^^^Hg results. Validity of the RCC 
results are countenanced by comparing our calculations at the zeroth order Dirac-Fock method and 
all order random-phase approximation with the other groups. In view of the current experimental 
progress in the ^^®Ra EDM measurement, our high accuracy calculations are very useful to yield 
limits on Ct and S in this atom. Thusly, it is possible to realize more accurate limits on the 
electron-quark T-PT interaction and Oqcd parameter of particle physics in future. 

PACS numbers: 


The experimental group at Argonne National Labora¬ 
tory (ANL) has just communicated the first measurement 
of permanent electric dipole moment (EDM) of ^^®Ra [l| 
after delineating its progress gradually in the last several 
years M- They obtain an upper limit on the atomic 
EDM (c^a) as |(iA(^^^Ra)| < 5.0 x 10“^^|e|cm (at 95% 
confidence). Though this limit is not competitive with 
the present current best limit available from ^®®Hg Q as 
|dA(^®®Hg)| < 3.1 X 10“^®|e|cm (at the same 95% con¬ 
fidence level), the ANL group portray their credence to 
surpass the limit on ^^(^^^Ra) over (iA(^®®Hg) sooner by 
making rapid improvement in their measurement @ . In¬ 
deed, viability of such claims are conceivable from both 
the theoretical and experimental prospectives. From the 
theoretical notion, the octupole deformation can enhance 
EDM observation in ^^®Ra by a factor larger than 10^ 
compared to ^®®Hg [3, Q . Other theoretical advantages 
that favor to carry out ^^®Ra EDM measurement are, 
it has larger nuclear charge (Z) compared to ^®®Hg and 
uncertainty due to the contribution from the octupole 
moment is still absent owing to its nuclear spin I = 1/2. 
On the experimental scope, cold-atom techniques to mea¬ 
sure Larmour precision have been developed for ^^^Ra 
atoms Si that are least sensitive to the systematics 
[ij . Moreover, the ANL research group also asserts to 
bring about significant improvement in the statistical un¬ 
certainty by increasing the number of atoms to 10® by 
producing them using the facility for rare isotope beam 
(FRIB) with the measurement time of 100 days Q. 

The EDM of a diamagnetic atom like ^^®Ra is sen¬ 
sitive to the parity (P) and time reversal (T) violating 
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(P,T-odd) electron-nucleon tensor-pseudotensor (T-PT) 
and nuclear Schiff moment (NSM) interactions. The 
NSM originates primarily due to the distorted charge 
distribution inside the nucleus caused by the P,T-odd 
interactions among the nucleons mediated by the neu¬ 
tral pion (7r®-meson). The other possible source to the 
^^®Ra EDM comes from the possible EDM or chromo- 
EDM of either up ((i„) or down [da) quark through the 
self-interactions in the 7r-meson loops [Ill) ■ Thus at 
the first hand if EDM of an atom is being observed, it 
will be a clean signature of P,T-odd interaction inside 
the atom [l^, . The second but the most vital conse¬ 

quence is, it will paramount to observing semi-leptonic 
CP-violation via the CPT theorem 15| and can explain 
certain new physics beyond the standard model (SM) of 
particle physics. In the SM, the complex (5-phase appear¬ 
ing in the Cabibbo-Kobayashi-Maskawa (CKM) mecha¬ 
nism is responsible for the hadronic CP-violation. Till 
date the CP-violation only in the hadronic sector has 
been observed [3) [l3| • The SM also predicts EDMs of 
elementary particles, but with sufficiently smaller values. 
However, it is largely believed that SM is an intermedi¬ 
ate manifestation of a complete theory. Such insistences 
are strongly endorsed by observation of finite masses 
of neutrinos, matter-antimatter asymmetry in the Uni¬ 
verse, existence of dark matter (indirectly) etc. 3, 1^ . 
Many models extension to SM have been propounded as 
attempts to perceive explanation to these unsettled in¬ 
quests. Some of the prominent models capable of explain¬ 
ing physics of the above problems are the left-right sym¬ 
metric, simersymmetric (SUSY), multi-Higgs etc. mod¬ 
els [III, [ifl [HI . These extensions to SM try to offer ad¬ 
ditional interactions, new particles and also incorporate 
other possible sources of CP-violation that can enhance 
EDMs of the elementary particles as well as of the com¬ 
posite systems up to the level as par with the current 
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sensitivity achieved by the atomic experiments. Thereby 
even if the limits of the T-PT coupling coefficients or 
NSM are improved from the atomic EDM studies, it can 
constrain to certain beyond SM models. In fact, NSM is 
related to Oqcd parameter of the particle physics |lll.l23| 
and improved value on 9qcd can annex to new sources of 
CP-violation apart from the (5-phase of SM which may be 
be sufficient enough to uphold the validity of those mod¬ 
els that are able to justify simultaneously large value of 
Oqcd and matter-antimatter asymmetry of the Universe. 
Therefore the ongoing efforts to improve limits on the ex¬ 
tracted parameters from the atomic EDMs, even with the 
null results, are still of pragmatic use. 

With the current development in the EDM measure¬ 
ment of ^^^Ra, it is now absolutely necessary to perform 
the required calculations accurately so that better lim¬ 
its on the above discussed parameters can be deduced. 
Recently, we have developed and employed a number of 
many-body methods in the relativistic framework to cal¬ 
culate these quantities more reliably among which the rel¬ 
ativistic coupled-cluster (RCC) methods are more pow¬ 
erful In this paper, we present very accurate cIa 

calculations for ^^®Ra due to both the T-PT and 

NSM interactions. We observe unusually large 

correlation effects in the determination of these proper¬ 
ties, but exhibiting strong cancellations for which the fi¬ 
nal results are reduced substantially compared to the val¬ 
ues reported previously using the methods having facile 
approximations. In order to demonstrate the same, we 
also compare our results at the intermediate calculations 
considering the Dirac-Fock (DF) method and random- 
phase approximation (RPA). We, then, appraise the role 
of the physical effects (specifically the non-RPA correla¬ 
tions) responsible for reducing the results that are usually 
appearing through the higher order perturbation theory. 
Moreover, we also perform calculation of the dipole po¬ 
larizability (ad) of ^^^Ra by taking the electric dipole 
(El) operator, which has identical angular momentum 
and parity selection rules with the T-PT and NSM in¬ 
teraction Hamiltonians, to compare its value with those 
of the previous calculations. This can provide additional 
insights into the potentials of various methods to yield 
accurate results. 

The P,T-odd electron-nucleus (e-N) interaction Hamil¬ 
tonians due to T-PT coupling and NSM are given by 
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respectively, with Gf is the Fermi coupling constant, Ct 
is the T-PT coupling constant, Sn= {Tn)j is the Pauli 
spinor of the nucleus, 7 d represents for the Dirac matri¬ 


ces, Pn{r) is the nuclear charge density, S = Sj is the 
NSM and B 4 = drr^pnir). 

The ground state wave function (|'I'q°^)) of ^^®Ra is 
evaluated using the relativistic coupled-cluster (RCC) 
method considering the dominant atomic Dirac-Coulomb 
(DC) Hamiltonian (Hoc) by expressing j'l'o) = 1‘ho) 

with the DF wave function |$o) and the RCC opera¬ 
tor that formulates excitation configurations acting 
upon |<I>o) to account for the correlation effects. For the 
determination of dA and aa, we also determine the first 
order wave functions (I'I'q^^)s) with respect to |'I'g°^) by 
considering the respective P,T-odd interaction Hamilto¬ 
nians and D operator. In the RCC framework, we express 


|vl/W) (3) 


where is also the RCC operator similar to but 
generating all possible odd parity excitations from |<i)o). 
For the calculations of both and we only 

allow all possible single and double excitation configura¬ 
tions (CCSD method). In this approach, both dA and ad 
(b^h a re symbolically denoted as X) is determined by 


A = 2($o|(i5^°^r(i)),,|$o) (4) 

where cc stands for the terms those are closed and con¬ 
nected and D = . In our previous work 

[ 2 ^, we have explained the procedure by which we take 
into account contributions from the non-truncative se¬ 
ries D = D 6 ^''°^ by an iterative procedure. We 

present later results even with the intermediate trunca¬ 
tion in this series as CCSD*^^^ method involving k num¬ 
ber of operators to highlight relevance of including 
higher order terms for accurate determination of dA and 
ad (we shall refer to our CCSD(°°^ method to simply as 
CCSD method afterwards). In fact, our CCSD method 
includes contributions from lower order many-body per¬ 
turbation theory (MBPT) and RPA. Therefore, the dif¬ 
ferences observed between the RPA and CCSD results 
are mainly because of the non-RPA contributions. Few 
important lower order MBPT diagrams representing non- 
RPA contributions are shown in Fig. [TJ 

We optimize size of the basis to be able to carry out 
CCSD calculations within the available computational 
resources by performing RPA calculations (that requires 
less computational time and memory compared to the 
CCSD method) considering a different combination of 
low-lying single particle orbitals of ^^^Ra from the s, p, 
d and g angular momentum symmetries. The RPA re¬ 
sults from few selective set of bases are given in Table 
HI and as seen, taking orbitals from beyond basis set IV 
do not change the results in the present interest of ac¬ 
curacy. Our calculations from different methods using 
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FIG. 1: Few important non-RPA diagrams from the MBPT 
method. Lines with up and down arrows are the occupied and 
unoccupied orbitals, respectively. Hi„t corresponds to either 
or D in the evaluation of and 

ad respectively. 


TABLE I: Demonstration of convergence in the q^, 
and results, in units of coq, (10“^*’C'T{o'n)|e|cm) and 

[S'/lei/m®]|e|cm) respectively with the finite bases, 
considering orbitals from the s, p, d, f and g angular mo¬ 
mentum symmetries using RPA. 


Basis 

S 

P 

d 

/ 

9 

OLd 

Results 

jTPT 

“A 

iNSM 

“A 

I 

10 

9 

8 

5 

4 

295.49 

-11.60 

-5.58 

II 

12 

11 

10 

7 

6 

296.79 

-14.81 

-7.20 

III 

14 

13 

12 

9 

7 

296.84 

-16.10 

-7.84 

IV 

17 

16 

15 

11 

9 

296.85 

-16.66 

-8.12 

V 

18 

17 

16 

12 

10 

296.85 

-16.66 

-8.12 


basis set IV are reported in Table HT] To estimate the un¬ 
certainties in our calculations, we estimate contributions 
from the orbitals of higher angular momentum symme¬ 
tries that were not considered in the calculations, fre¬ 
quency independent Breit interaction and leading order 
quantum electrodynamics (QED) corrections (through 
the Uehling, Wichmann-Kroll and self-energy accounting 
potentials as described in [^) using RPA. These contri¬ 
butions are also quoted separately in Table HIl which show 
that the neglected contributions may not be larger than 
1 %. 

Though we have estimated uncertainties by taking into 
account various truncations in the methodologies of our 
numerical calculations, however it would at least be de¬ 
sirable to compare calculations of a physical quantity us¬ 
ing the employed methods with its experimental result to 
test their capabilities for producing accurate result. As 
said before, evaluation of both dA and aa have some simi¬ 
larities in the sense that their determination entails same 
angular momentum and parity selection criteria. On that 
account, calculations of dA and ad are often performed 


together 13-25, 28- 


albeit it demands for different 
radial behavior of the wave functions for their estima¬ 
tion. Yet the experimental result of ad in ^^®Ra to be 
realized, but many exclusive calculations using variants 
of RCC method report its value [lll - fsljl . Borschevsky et 
al. have used a CCSD method including important par¬ 
tial triples excitations perturbatively (CCSD(T) method) 
with the DC Hamiltonian and Gaunt term of the Breit in¬ 
teraction [Hj . The fundamental difference between their 
calculation with ours is that Borschevsky et al. use a 
finite field approach through a molecular code while we 
use the expectation value determination method describ¬ 
ing wave functions in the spherical coordinate system. 
Consequently, Borschevsky et al. obtain large DF value 
compared to their CCSD(T) result. On the other hand, 
our CCSD result is larger than the DF value. Neverthe¬ 
less, the final results in both the cases (that are suppose 
to converge towards the exact value) seem to be match¬ 
ing with the difference of about 2%. On the other hand 
there are also two more calculations available using the 
methods similar to Borschevsky et al. but considering 
the scalar relativistic Douglas-Kroll (DK) Hamiltonian 
and small-core ten-valence electron scalar-relativistic 
pseudopotentials (ARPP) [s^l. Very good agreement be¬ 
tween our CCSD result with the CCSD(T) values for ad 
of ^^^Ra entails that our method, which is really an ardu¬ 
ous procedure to implement in the spherical coordinate 
system, is adequate for the accurate determination of dA 
and ad. 

After having a satisfactory comparison between the ad 
values from the CCSD methods, we also analyze both the 
dA and ad results from Table El among other approxi¬ 
mated methods using which dA of ^^®Ra were previously 
evaluated. Among the other previous theoretical studies, 
the combined conhguration interaction and leading order 
many-body perturbation theory (CI-I-MBPT method) 
with some contributions through time-dependent DF 
method (equivalent to our RPA) employed by Dzuba and 
Coworkers seem to include more phy sical effects in the 
evaluation of both dA and ad [23, l30| ■ In this approach, 
they have used a V^~‘^ potential, N = 88 is the to¬ 
tal number of electrons of ^^®Ra, to generate the single 
particle orbitals in contrast to our potential. Un¬ 
like our methods, the correlation effects from different 
electrons are not accounted in equal footing in these cal¬ 
culations. Also, Latha and Amjith have employed RPA 
with potential to report dA due to T-PT interaction 


and ad [29|. Results at RPA level match well between 
all these calculations. The CI-fMBPT results of Dzuba 
and Coworkers, however, give larger dA but smaller ad 
values compared to the RPA calculations. On the oth- 
erhand, correlation effects through CCSD method (espe¬ 
cially the all order non-RPA contributions) reduce signif¬ 
icantly both dA and ad with respect to the RPA results. 
It should be reminded that our CCSD method already 
includes RPA result. The and d^^^ values from 
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TABLE II: Estimated ad, and (in same unit as of 

Tablelil) results for ^^®Ra from various calculations. CCSD^°°^ 
values are considered as our final results and contributions 
from the Breit interaction {5b), QED correction (Sqed) and 
truncated basis (Sbasis) are accounted for estimating uncer¬ 
tainties to the CCSD results. 


Method of 

This work 


Others 


Evaluation 

ad 

iTPT 

dA 

iNSM 

dA 

(^d 

jTPT 

dA 

jNSM 

dA 

DF 

204.13 

-3.46 

-1.85 

204.2 

200.6*’ 

293.4'’ 

300.57'* 

299.52" 

-3.5“ 

-1.8“ 

RPA 

296.85 

-16.66 

-8.12 

291.4*' 

-17“ 

-16.59*’ 

1 

00 

CO 

p 





297.0^ 


-8.5^ 

CI-bMBPT 




229.9-^ 

-18“ 

-8.8“ 

CCSD^ 

253.04 

-10.40 

-6.94 




CCSD^ 

242.02 

-9.49 

-6.52 




ccsd(°°) 

CCSD(T) 

247.76 

-10.04 

-6.79 

242.8" 

248.56'* 

251.12" 



Sb 

0.19 

0.06 

0.06 




Sqed 

-0.43 

-0.16 

-0.07 




^basis 

-0.03 

-0.08 

-0.05 





“[ 2 ^ 

*’[2^ 

Corrections from the Gaunt term incorporated. 
DK operator is used in the calculations. 

ARPP is considered in the calculations. 

^[30] Corrections from the RPA included. 


the CCSD method are reduced by 45% and 23%, respec¬ 
tively, compared to the CI-I-MBPT results for ^^®Ra even 
though difference between the ad value from both the cal¬ 
culations are only about 10%. We also present CCSD^^^ 
and CCSD^"^) results in Table HIl to demonstrate the role 
of the correlation effects to improve the results gradually. 

Combining our final CCSD results with the mea¬ 
sured dAi^^^Ra) Q, we get an upper bound on NSM 
as S'<7.4 X 10“®|e|fm^. Similarly with the knowledge of 
(an) from the nuclear calculation, upper bound on Ct 
can also be obtained. Two sophisticated nuclear struc¬ 
ture calculations by considering the octupole deformed 
Wood-Saxon^tential [s^ and the odd-A Skyrme mean 
field theory [33 have been carried out to describe the 
P,T-odd interactions in the ^^^Ra nucleus. These cal¬ 
culations parameterize S in terms of the pion-nucleon- 
nucleon (ttNN) couplings. However, they provide dif¬ 
ferent combinations of the ttNN couplings to express S. 
Recently Engel et al. in their recent review give the best 


value for S from these two calculations in terms of the 
ttNN couplings as [1^ 

S = 13.5[-1.55o + b.Ogi - 4 . 0 ^ 2 ] |e| W, (5) 


where 5 i=o,i ,2 are the isospin components of the P,T-odd 
ttNN coupling constants. Thus, combining the above 
result with our limit on S, we infer bounds on these 
couplings as |go|<3.6 x 10“^ and |gi|<9.1 x 10“®. Fur¬ 
ther from the relations |go| = O.O18(7)0 qcd l37| and 
|gi| = 2 X —Jd) (Hi) we extract the upper limits 

as |0 qcd|<2.O X 10“® and |J„ — dd\ < 4.6 x 10“^^|e|cm. 
Though these bounds are not competitive at present 
against those obtained using the EDM study of ^®®Hg 
[23j, however as said before rapid improvement in the 
EDM measurement of ^^®Ra is anticipated [ij. It, there¬ 
fore, appears very likely that limits on the above quan¬ 
tities can become more stringent with the upcoming im¬ 
proved EDM measurement of ^^^Ra and when the nuclear 
calculations are advanced further. These limits would 
definitely shed light on the new physics beyond the SM. 

We thank Professor B P Das for encouraging us to 
carry out this study. We are grateful to Professor Z -T 
Lu for informing us their ^^®Ra EDM experimental result. 
The computations were carried out using 3TFLOP HPC 
cluster of PRL, Ahmedabad. 
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